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SHORT  WAVELENGTH  STABILIZATION  OF  THE  GRADIENT  DRIFT  INSTABILITY 

DUE  TO  VELOCITY  SHEAR 


I.  INTRODUCTION 

The  gradient  drift  Instability  (also  known  as  the  E  x  B  Instability  or 

the  E  x  B  gradient  drift  Instability)  has  been  of  Interest  to  space  and 

plasma  physicists  for  two  decades.  The  Instability  was  Initially  studied 

by  Simon  (1963)  and  Hoh  (1963)  to  explain  turbulence  observed  In  laboratory 

experiments.  Since  then  the  dominant  Interest  In  the  Instability  has  been 

In  regard  to  natural  Ionospheric  Irregularities  [e.g.,  equatorial 

electrojet  Irregularities  (Farley,  1979;  Fejer  and  Kelley,  1980)  and  high 

latitude  irregularities  (Keskinen  and  Ossakow,  1982)]  and  to  man-made 

Ionospheric  Irregularities  [e.g.,  plasma  cloud  releases  (Llnson  and 

Workman,  1970)].  The  instability  can  arise  In  a  weakly  colllslonal  plasma 

which  contains  a  density  gradient  that  Is  orthogonal  to  the  current  flow 

generated  by  an  ambient  electric  field  perpendicular  to  the  magnetic 

field.  Depending  upon  the  ratios  v  /ft  and  v./ft.,  where  v  ....  Is  the 

e  e  l  l  e(l) 

electron  (ion)-neutral  collision  frequency  and  ft  .  Is  the  electron  (Ion) 

ei  t ; 

cyclotron  frequency,  an  ambient  electric  field  can  produce  two  types  of 

current  in  a  weakly  colllslonal  plasma.  For  the  case  of  v  /.ft  «  l  and 

e  e 

Ui^i  ^  a  Pedersen  current  Is  produced  by  the  Ions  in  the  direction  of 
the  electric  field;  when  v  /ft  «  1  and  v./ft.  »  1,  a  Hall  current  Is 
produced  by  the  electrons  In  the  direction  orthogonal  to  the  electric 

field.  The  Pedersen  current  driven  instability  is  relevant  to  F  region 
Irregularities,  while  the  Hall  current  driven  instability  is  relevant  to  F, 
region  irregularities.  This  letter  will  discuss  the  influence  of  velocity 

shear  on  the  gradient  drift  instability  driven  by  the  Hall  current. 

Manuscript  approved  February  9,  1983. 
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The  motivation  for  this  paper  is  two-fold.  First,  studies  of  the 
Rayletgh-Taylor  instability  (Guzdar  et  al.t  1982)  and  the  Pedersen  current 
driven  gradient  drift  (or  E  x  B)  instability  (Perkins  and  Doles,  1975;  Huba 
et  at. ,  1983)  indicate  that  velocity  shear  can  have  an  important  effect  on 
these  instabilities.  Namely,  the  stabilization  (or  suppression)  of  the 
short  wavelength  modes  which  causes  a  preferential  excitation  of  a  longer 
wavelength  mode.  Thus,  it  is  anticipated  that  velocity  shear  will  have  a 
similar  effect  on  the  Hall  current  driven  gradient  drift  instability,  and, 
in  fact,  we  show  that  this  is  the  case.  Second,  recent  studies  of  long 
wavelength  turbulence  in  the  equatorial  electro jet  (Kurieki  et  al.,  1982; 
Pfaff  et  al.,  1982)  indicate  that  a  long  wavelength  mode  with  X  few  km  is 
dominant  during  type  T  conditions.  Thus,  inclusion  of  velocity  shear  in 
the  plasma  configuration  of  the  electrojet  may  account  for  this 
observation. 
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II.  THEORY 


The  plasma  configuration  and  slab  geometry  used  In  the  analysis  are 
shown  In  Fig.  1.  The  ambient  magnetic  and  electric  fields  are  In  the  z  and 

A  A 

-x  directions,  respectively  (l.e.,  §  a  ez  an<*  5  “  -  Eq(x)  ex^'  ^he 

density  and  electric  field  are  assumed  to  be  Inhomogeneous  In  the  x 

direction.  The  collision  frequencies  are  such  that  v  «  2  and  v.  »  Q. 

e  e  I  l 

where  Is  the  a  specie-neutral  collision  frequency  (assumed  to  be 

constant)  and  ft  *  |e  |B„/ra  c  is  the  cyclotron  frequency  of 
a  a  0  a 

the  a  species.  The  electrons  can  be  considered  to  be  strongly  magnetized, 
while  the  Ions  are  unmagnetized.  Thus,  the  electrons  experience  a 

A 

sheared  E  x  B  drift  In  the  y  direction  (V.(x)  -  -  cE_(x)/B_  e  )  while  the 
~  ~  vo  U  0  y 

ions  do  not.  This  differential  motion  between  the  electrons  and  ions 
produces  a  Hall  current. 

We  assume  that  perturbed  quantities  vary  as  6p  ~  5p(x)exp[ l(k^y  -  at)] 
where  ky  is  the  wave  number  In  the  y  direction  and  u>  ■  u>  +  iy  so 
that  y  >  0  implies  wave  growth.  We  neglect  perturbations  along  the 
magnetic  field  (k  •  B~  ■  0)  so  that  only  the  two  dimensional  mode  structure 

v  vU 

in  the  xy  plane  Is  obtained.  We  assume  so  that  electron  inertia 

terms  can  be  neglected.  An  important  feature  of  our  analysis  Is  that  a 
nonlocal  theory  Is  developed.  That  Is,  the  mode  structure  of  the  potential 
In  the  x  direction  (the  direction  In  which  the  density  and  electron  drift 
velocity  are  assumed  to  vary)  is  determined  by  a  differential  equation 
rather  than  an  algebraic  equation  obtained  by  Fourier  analysis.  This  Is 
crucial  to  the  analysis  since  previous  studies  of  interchange  instabilities 
using  a  sheared  velocity  flow  equilibrium  [Perkins  and  Doles,  1975;  Guzdar 
et  al.,  1981;  Huba  et  al.,  1983]  have  shown  that  a  nonlocal  theory  is 
necessary  to  demonstrate  the  stabilizing  Influence  of  velocity  shear. 
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The  fundamental  equations  used  in  the  analysis  are  continuity  and 
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momentum  transfer  In  the  frame  of  reference  of  the  neutrals; 


3n 


~  +  V  .  (n  V  )  =■  0 
3t  v  ct'-cr 


1  1 

-  £_  (e  +  -  V  x  - —  7n  -  v  V 

m  c  ~e  -0-'  nm  e~e 


dV.  T. 

v  j  g  i 

=  — —  E  -  -±-  Vn  -  \> .  V 
dt  m^  x  nm^  l~l 


(1) 

(2) 

(3) 


where  a  denotes  partLcle  species  (e:  electrons;  l:  Ions)  and  other 
variables  have  their  usual  meanings.  From  Eqs.  (2)  and  (3)  we  find  that 
the  equilibrium  drifts  of  the  electrons  and  tons,  to  first  order 
In  ve/fte  and  are  given  by 


V  ■  7T-  (vj  “  V.  )e  +  (V.  -  V.  )e 
-e  n  ^  d  de-'  x  v  d  de'  \ 


and 


V.  =■  -  —  (V.  -  V ,  )e 
xi  ^  d  di'  } 


(4) 


(5) 


where  *  -  cEq/Bq  Is  the  E  x  B  velocity  and  =■  (cT^/e^B^DJln  n/3x  is 
the  diamagnetic  drift  velocity  of  species  a.  An  equilibrium  relationship 
between  n  and  can  be  obtained  by  substituting  Eqs.  (4)  and  (5)  into  Eq. 
(1)  and  setting  3n^/3t  *  0.  Assuming  »  V^,  the  equilibrium  relation 
between  n  and  is 
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n(x)V,(x)  =  constant 

U 


(6) 
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I 

i 
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where  we  have  neglected  the  x  variation  of  the  collision  frequencies.  We 
hasten  to  add  that  Eq.  (6),  although  a  self-consistent  relationship  within 
the  context  of  Eqs.  (1)  -  (3),  does  not  necessarily  model  the  equatorial 
electrojet  or  auroral  plasmas. 

We  linearize  Eqs.  ( 1 )— ( 3 )  by  assuming  n  -  n  +  fin  ,  V  =■  V  +  fiV 

CL  G  01  G  Ot 

and  E  ■  Eq  -  V,J>  where  $  is  the  perturbed  electrostatic  potential.  Since  we 

assume  v  «  Q  ,  v,  »  3,,  and  V.  »  V.  ,  we  take  the  equilibrium  drifts  to 
e  e  i  l  d  dot 

be,  to  lowest  order, 


V 

~e 


V  (x)e 
d  y 


-  cE0(x)/Bq  ey 


(7) 


V 


i 


0. 


(8) 


I 

I 

I 

i 

►n 

u 


The  perturbed  velocities  are  given  by 


Ik 


oV 


[ - L  (®_  6$  -  v2  — -1  +  -S.  f  ~  4,  -  V2  — -)]e 

l  '■m  v  e  n  '  Q  Q  vmy  enJJx 


fin 


v  ik 
e  k  re 


5n 


e  e 


ik  fin  v  ik  5n 

+  [-i  (i-  6q>  -  v2  — — )  +  -—2-  (  —  <}>  -  v2  -£■)]■ 

1 3  '-m  y  en-'  3  3  ^  ra  y  en  J 

q  G  C  6  6 


ik  fin, 

—  (—  <|>  +  v2  - •) 


(9) 


(10) 


1/2  * 
where  v  *  f T  /m  ]  Is  the  thermal  speed  of  species  a  and  k  *  k  e  + 
a  y  a  aJ  ~xx 

a  a  a 

kyey  »  -  i  3/3x  ex  +  k yoy,  i»e.,  ^  is  an  operator.  In  writing  Eqs.  (9) 

and  (10)  we  have  assumed  k2L2  »  1  and  k2L2  »  1  where  L  ■  (3  £n  n/3x)  \ 

y  n  x  n  n 

so  that  we  are  only  considering  unstable  modes  which  have  wavelengths  small 
compared  to  the  scale  length  of  the  density  Inhomogeneity. 

We  substitute  Eqs.  (9)  and  (10)  into  Eq.  (1)  and  obtain 
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f-  if  to  -  k  V.)  +  k2p2v  +  ik  p  v  /L  1  —S. 
L  ^  y  dJ  e  e  y  e  e  nJ  n 


f  lk  )-  +  k2  -£■]  —  i-  4 

1  1  S,  V  ■«  \ 


Oil  j 

[“(“  +  ivi>  "  k2vil  —  =  ^  k2*’ 


where  ■  v>^e  an<*  =  ^  n/<*x)  •  Assuming  quasi-neutrality 

(Sn  *  5n. )  we  find  that  Eqs.  (11)  and  (12)  reduce  to 
s  1 

v  v  k  ik 

L  clk4^  +  [  (to  -  k  1^)8^  "iw(w  +  ivt)  ~  -  c  g  ]  k2  4>  +  =  0 


finally,  we  assume  k2  «  k2  and  expand  k4  about  the  small  parameter 

x  y 

k2/k2,  retaining  terms  to  second  order  in  kx.  We  then  make  the  identifi- 
x  y 

cation  k2  =  -  3 2 / 3 x2  and  arrive  at  the  mode  equation 


-  k2  Q(x)*  -  0 
3x2  y 


where  Q(x)  ■  q(x)/p(x)  and 


p(x)  =  to  -  kyvd(x)  -  kypi[  ^kyLn'  ”  2l^e]Ql 


q(x)  -  to  -  k  Vd(x)  +  [to(a>/ni  +  ivt)  -  kyp^nj^k  L  )“  -  ivg 


where  p.  -  c  Al,,  v.  -  \>./n,  and  v  -  v  /3  . 
l  si  t  tl  e  ee 
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III.  RESULTS 


In  order  to  understand  the  Influence  of  velocity  shear  on  the  gradient 
drift  instability  we  contrast  two  cases.  In  the  first  case  we  assume  Eq  * 
constant  (i.e.,  Vj  =  constant;  no  shear)  and  in  the  second  case  we  consider 
Eq  *  constant  (i.e.,  =»  V^x);  shear). 

A.  No  shear  (Vj  =  constant) 

We  assume  Eq  =  constant  so  that  =  constant,  and  use  local  theor> 

previous  authors  have  considered  (e.g.,  Roglster  and  D^Angelo,  1970;  _Si 

et  al. ,  1973),  i.e.,  Eq.  (6)  is  violated.  The  local  approximation  assumes 

k2L2  »  1  and  k2L2  »  1,  i.e.,  the  wavelengths  of  the  unstable  modes  are 
x  n  y  n  ° 

small  compared  to  the  density  inhoraogeneity .  The  dispersion  equation  (Eq. 
(13))  can  be  written  as 


“( 1  +  Vi)  -  Vd  "  /  ET  +  57  (“2  '  k2c|)l  •  0 


(17) 


where  we  have  assumed  vp  »  1/kL^  for  simplicity.  Assuming  |yl  «  U  | ,  we 
find  the  solution  of  Eq.  (17)  to  be 


Y 


o  v, 

e  1 


•a  k 


JL  JL  JL 


.  -  -  L  -  k  kL 

1  +  v  v.  v  n 

e  i  e 


Vi 


(18) 


(19) 


Equations  (18)  and  (19)  agree  with  earlier  results  (Rogister  and  D'Angelo, 
1970;  Sudan  et  al.,  1973)  and  describe  both  the  Farley-Buneman  and  gradient 
drift  instabilities.  The  Farley-Buneman  Instability  can  occur  when 
u)  >  kcg  or  Vj  >  c8(l  +  vive)*  0n  the  other  hand,  the  gradient  drift 
instability  is  excited  when 


I 


s 


i 


* 

i 


i 


* 


k  (cEn/B)  .  k  . 
yv  0  1  y  l 

^  v  ic  kL 

1  +  v .  v  v  n 

i  e  e 


>  0 


(20) 


for  kp ,  «  1.  When  V.  <  c0(l  +  v.v  )  and  kp,  -  1,  the  gradient  drift  modes 
i  as  i  e  i 

are  stabilized  because  of  electron-neutral  collision  (l.e.,  diffusion 
damping).  The  important  condition  represented  by  Eq.  (20)  Is  that  the 
electric  field  and  density  gradient  be  in  the  same  direction  for 
instability. 

We  plot  the  growth  rate  y  /ft.  vs.  k  p  (and  kvL)  in  Fig.  2  (curves  A 

q  i  y  I  y 

and  B)  for  the  case  of  no  shear.  Equation  (13)  is  solved  numerically  for 

the  following  parameters:  Vj/cs  =  "1*5  (curve  A)  and  V^/Cg  =  -0.2 

(curve  3),  p^/L  =  10  \  Ln/L  =  1.0,  =  10  t  and  ~  25.0. 

Here,  is  the  growth  rate  maximized  with  respect  to  kx/k y.  We  note  the 

following.  First,  Ym  *  k^  for  k^p^  $  0.04  (or  kyL  <  40)  in  both  curves  A 

and  B.  Second,  y  ~  constant  for  0.04  <  k  p,  <  0.40  (or  40  <  k  L  <  400). 

And  finally,  y^  increases  sharply  for  k^p^  >  1  when  Vj/cs  =  -1.3  which  is 

due  to  the  Farley-Buneraan  instability  (curve  A);  y^  becomes 

negative  (y  <  0)  for  k  p,  -  0.65  when  Vj/cc  =  -0.2  which  is  due  to  the 
in  y  1  n  ?> 

colltsional  electron  damping  of  the  gradient  drift  Instability  (curve  3). 

B.  Shear  (V,  t  constant) 
d 

We  assume  F,q  a  constant  so  that  =  V^(x)  and  a  sheared  velocity  flow 
exists  In  the  plasma.  As  in  the  case  of  the  E  x  B  Instability  [Perkins  and 
Doles ,  1975;  Huba  et  al.  ,  1983]  local  theory  Is  inadequate  to  properly 

describe  the  linear  properties  of  the  mode.  We  solve  Eq.  (14)  numerically 
using  a  finite  difference  scheme  [i.e.,  the  Nuraerov  method  (Gladd  and 
Horton,  1973)]  to  obtain  eigenvalues  and  eigenfunctions.  The  boundary 
conditions  used  are 


i 
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Plot  of  the  maximum  growth  rate  vs.  wavenvimber 

[y  /Q,  vs.  k  L  (k  p,)l  for  the  cases  of  no  shear  (A  and  B)  and 
ml  y  jr  i  1 

shear  (A'  and  B');  and  of  the  phase  velocity  vs.  wavenumber 

[V  .  *  w/k  c  vs.  k  L  (k  p,)j,  denoted  by  circles,  for  the 
ph  y  s  y  y  1  '  J 

sheared  case.  See  the  text  for  a  description  of  the  parameters. 


with  the  sign  of  the  WKB  solution  cho  .-.n  such  that  a  damped  solution  Is 
obtained  in  the  limit  |x|  + 

The  density  profile  used  In  the  analysis  Is  given  by 


1  +  e  tanh(x/L) 
1  -  e 


(22) 


where  0  <  e  <  l,  L  characterizes  the  width  of  the  boundary  layer,  and 
nQ  -  n(x  »  -  ®).  We  assume  e  *  0.8  for  the  results  presented  so  that 


max 


r  3  in  n' 

^  3x  -*raax 


x/L 


-0.55 


(23) 


The  equilibrium  velocity  profile  chosen  to  satisfy  Eq.  (6),  l.e., 


Vd(x)  "  Vd0n0/n(x)  (24) 

where  Vd0  -  Vd(x  »  -  «). 

In  Fig.  2  we  plot  Ym/0^  vs.  kyL  (or  k^p^)  based  upon  Eqs.  (14),  (22) 

and  (23)  for  the  following  parameters:  VdQ/cg  “  “7.5  (curve  A')  and 

Vd0^cs  “  (curve  B'),  p^L  *  10  3,  *  10  2,  and  =  25.0. 

Here,  is  the  growth  rate  maximized  with  respect  to  mode  number,  which  Is 

analogous  to  kx.  We  point  out  that  the  mode  structure  of  the  unstable 

waves  (l.e.,  the  eigenfunction)  Is  not  localized  about  the  position  of 

maximum  density  gradient  (x  ■  -0.55  L)  but  rather  at  x  =  0.  The  local 

drift  velocity  at  x  ■  0  is  approximately  V,„  =  Vjrv/5  so  that  V,„  =  -1.5 

dJt  au  dx 


(curve  A')  and  Vd  =  -0.2  (curve  B') 


Thus,  the  local  plasma  conditions 


are  roughly  the  same  between  curves  A  and  A',  and  between  B  and  B'.  The 

only  difference  being  that  the  primed  curves  contain  the  effects  of 

-2 

velocity  shear.  In  the  region  k^L  <  10  (or  kyP^  <10  )  it  Is  clear  that 

velocity  shear  has  no  effect  on  the  gradient  drift  Instability  since  curves 

-2 

A  and  A',  and  B  and  B'  coincide.  However,  for  kyL  >  10  (or  k^p^  >10  ) 

the  primed  and  unprlmed  curves  diverge.  We  find  that  the  gradient  drift 
Instability  achieves  maximum  growth  at  kyL  =  20  and  is  stable  (y  <  0)  at 
kyL  »  35.  The  value  of  kyL  for  which  the  mode  is  "shear  stabilized"  is 
more  than  an  order  of  magnitude  smaller  than  the  diffusive  cutoff 
associated  with  curve  B  (kyL  ~  650).  Thus,  velocity  shear  has  the  effect 
of  stabilizing  the  short  wavelength  gradient  drift  modes.  Note  that  the 
position  of  maximum  growth  is  independent  of  the  drift  velocity  V^.  Thus, 
the  wavelength  of  the  most  unstable  modes  is  dictated  by  L,  not  V^. 

We  also  plot  the  phase  velocity  of  the  waves  as  a  function  of  ky  in 
Fig.  2.  Specifically,  we  show  Vp^  *  u)/kycg  vs.  kyL  for  the  sheared  drift 
cases  (A'  and  B');  we  denote  Vp^  by  small  circles  in  Fig.  2.  For  kyL  <  20 
we  find  that  Vp^  «  ky,  while  for  kyL  >  20  that  ~  constant.  Unlike  the 
position  of  maximum  growth,  the  phase  velocity  of  the  waves  is  a  function 
of  the  drift  velocity  V(jj  the  larger  the  drift  velocity,  the  larger  the 
phase  velocity  of  the  waves.  For  the  parameters  considered  in  Fig.  2  we 
note  that  u>/k_c_  ~  0.1  for  V^/Cg  ■  -1.0  (B')  and  oj/k  c_  =  1.0  for 


y  s 

Wcs  “  "7*5  <A'>* 


y  s 
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IV.  DISCUSSION 


We  have  presented  a  nonlocal  analysis  of  the  gradient  drift 

instability.  The  new  effect  included  in  this  theory  is  the  allowance  for 

an  inhomogenous  electric  field  which  produces  a  sheared  drift  velocity, 

l.e.,  Vj  *  VjCx).  The  major  result  of  this  work  is  that  velocity  shear  can 

stabilize  the  short  wavelength  modes  of  the  Instability,  and  preferentially 

excite  a  longer  wavelength  mode  than  would  be  expected  from  local  theory. 

This  result  is  similar  to  the  short  wavelength  stabilization  of 

the  K  x  B  instability  [Perkins  and  Doles,  1975;  Huba  et  al.,  1983).  On  the 

other  hand,  the  long  wavelength  modes  appear  to  be  unaffected  by  velocity 

shear.  We  emphasize  that  we  have  not  explored  the  influence  of  velocity 

shear  on  the  Farley-Buneraan  Instability;  it  is  possible  that  y  in  curve  A' 

m 

of  Fig.  2  may  become  positive  for  k  p,  >  1.  We  also  mention  that  we  have 

y  i  - 

only  considered  linear  theory  and  that  a  nonlinear  cascade  of  wave  energy 
(Kesklnen  et  al.,  1979)  could  produce  wave  turbulence  in  the  regime 
kyP^  >  20  even  for  the  velocity  sheared  cases  considered  in  Fig.  2. 

We  discuss  the  possible  significance  of  these  results  to  the 
equatorial  and  auroral  E  regions.  Recently,  observations  of  long 
wavelength  irregularities  in  the  equatorial  F  region  have  been  made  (Kudekl 
et  al.,  1982;  Pfaff  et  al.,  1982).  These  long  wavelength  irregularities 
are  characterized  by  having  (1)  wavelengths  such  that  kL  -  23  where  L  -  3- 
30  km  (l.e.,  X  -  few  kms);  (2)  wavelengths  determined  by  the  gradient  scale 
length,  not  the  drift  velocity;  (3)  phase  velocities  Vp^  -  100  m/sec  < 
cs  ~  360  a/sec;  and  (4)  large  amplitudes.  Our  theory  is  consistent  with 
the  first  three  observations,  but  cannot  address  the  fourth  since  it  is 
only  a  linear  theory.  We  find  that  velocity  shear  stabilizes  the  short 
wavelength  modes  of  the  gradient  drift  instability  and  preferentially 


excites  a  mode  such  that  kyL  -  20;  that  the  scaling  of  the  dominant 

wavelength  Is  related  to  the  scale  size  of  the  system  and  not  to  the  drift 
velocity  (l.e.,  kyL  -  20  for  both  curves  A'  and  B');  and  that  the  phase 
velocity  of  the  dominant  wave  Is  such  that  Vp^  ~  w/k^  <  c8  for  Vj  <  cs. 

These  findings  suggest  that  velocity  shear  may  account  for  the  dominant 

long  wavelength  Irregularities  observed  In  the  equatorial  electro jet.  This 
conclusion  must  be  regarded  as  tentative,  though,  since  the  density  and 
velocity  profiles  used  In  this  analysis  [Eqs.  (22)  and  (24)],  based  upon 
Eq.  (6),  do  not  accurately  model  the  equatorial  electrojet.  However,  based 
upon  previous  work  (Guzdar  et  al.,  1982;  Huba  et  al.,  1983)  and  analysis  of 
Eq.  (14)  using  different  plasma  profiles,  we  find  that  the  fundamental 
conclusion  of  this  paper,  viz.,  velocity  shear  can  stabilize  short 
wavelength  modes  of  the  gradient  drift  Instability,  Is  not  sensitive  to  the 
plasma  profiles  chosen.  Nonetheless,  in  order  to  test  our  hypothesis 
regarding  the  equatorial  electrojet,  an  appropriate  plasma  equilibrium  Is 
needed  which  could  perhaps  be  supplied  by  experimentalists.  Also,  the 

actual  parameters  used  In  this  analysis  <«•*•»  Wcs  “  -  7*  5)  are  not 
appropriate  to  electrojet  conditions  but  are  chosen  to  contrast  the  role  of 
shear  to  the  role  of  drift  velocity. 

We  suggest  that  such  long  wavelength  irregularities  may  exist  In  the 
auroral  electrojet  and  could  be  responsible  for  the  long  wavelength 
structuring  of  auroral  Irregularities.  Discrete  auroras  are  observed  to 
have  two  distinct  size  scales:  (1)  one  Is  the  tens  of  kilometer  width 

associated  with  inverted  V  precipitation,  and  (ii)  the  other  is  the 
kilometer  widths  of  auroral  arc  elements  which  appear  to  be  imbedded  In  the 
inverted  V  precipitation  region  (e.g.,  Davis ,  1973).  The  formation  of 

inverted  V  arcs  can  be  understood  in  terms  of  the  magnetosphere-ionosphere 
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coupling  (e.g.,  Kan  and  Lee t  1980).  We  suggest  here  that  the  gradient 
drift  instability  discussed  in  this  paper  can  lead  to  the  development  of 
the  small-scale  arc  elements  within  the  inverted  V  precipitation  region. 
During  substorms,  the  ionospheric  plasma  density  in  the  auroral  arc  region 
is  greatly  enhanced  due  to  the  ionization  of  the  atmosphere  caused  by  the 
precipitation  of  energetic  electrons.  In  addition,  paired  electric  fields 
directed  toward  each  other  are  present  in  the  inverted  V  region. 
Therefore,  a  large  plasma  density  gradient  and  an  electric  field  gradient 
may  exist  in  the  arc  region,  which  will  result  in  the  long  wavelength 
gradient  drift  instability  as  obtained  here.  Take  an  inverted  V  region  to 
have  a  width  L  *  30  km.  The  eigenfunction  4>(x)  in  (14)  for  the  most 
unstable  mode  (k^L  =  20)  is  found  to  have  4-6  peaks.  Thus,  it  is  plausible 
that  the  large-scale  arc  will  be  filamented  into  4-6  arc  elements  as 
usually  observed  (Davis,  1978). 
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01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01 CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  HA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  CON  FOR  J.  MALIK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZINN 
OICY  ATTN  DOC  CON  FOR  P.  KEATON 
OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SAHDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 
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OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

(ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

01 CY  ATTN  R.  MOORE 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
OICY  ATTN  D.  CROMBIE 
OICY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  A  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICY  ATTN  R.  CRUBB 
OICY  ATTN  AERONOMY  LAB  G.  REID 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 


COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  CARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CALIFORNIA  AT  SAN  DISCO,  UNIV  OF 

P.O.  BOX  6049 

SAN  DIEGO,  CA  92106 


EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORCE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  ATTN  G.  MILLMAN 
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GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
OICY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTEKN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  6  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 


LAURAL, 

MD 

20810 

OICY 

ATTN 

DOCUMENT  LIBRARIAN 

OICY 

ATTN 

THOMAS  POTEMRA 

OICY 

ATTN 

JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS I AC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO. ,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  GLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  H  .RRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 
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MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  D.  SAPPENFIELD 
01CY  ATTN  R.  BOCUSCH 
01 CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONGMIRE 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  !1A  02173 

01CY  ATTN  IRVING  L.  KOFSKY 


R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01CY  ATTN  FORREST  CILMORE 
01CY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OiCY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIF.MAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
01CY  ATTN  RAY  L.  LEADABRAND 
01 CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  J.  PETERSON 
01CY  ATTN  R.  LIVINGSTON 
01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
01CY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OiCY  ATTN  W.P.  BOQUIST 

TOYON 

34  WALNUT  LAND 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 
OICY  ATTN  JOEL  GARBARINO 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
P1CY  ATTM  D/  Stockwell 
SNTF/1575 


VISIDYXE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

Dr.  S.  Ossakow  -  Code  4700  (26  copies) 
Code  4701 

Code  4780  (100  copies) 

DR.  P.  MANGE  -  CODE  4101 
DR.  E.  SZUSZCZEWICZ  -  CODE  4108 
DR.  J.  GOODMAN  -  CODE  4180 
DR.  P.  RODRIGUEZ  -  CODE  4108 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/ NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  R.F.  BENSON 
DR.  K.  MAEDA 
Dr.  S.  CURTIS 
Dr.  M.  DUBIN 

DR.  N.  MAYNARD  -  CODE  696 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  .ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  C.  JOINER 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMEKL 

GEOPHYSIICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.C.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  BARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 


UNIVERSITY  OF  CALIFORNIA, 
LOS  ANGELES 
405  HILLCARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE 'PARK,  KD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  CREENWALD 
DR.  C.  MENG 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 

UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  SPACE  SCIENCES 
P.O.  BOX  688 

RICHARDSON,  TEXAS  75080 
DR.  R.  HEELIS 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 


